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ABSTRACT 

We report the first detection of Blue Straggler Stars (BSS) in the bulge of the 
Milky Way galaxy. Proper motions from extensive space-based observations along a 
single sight-line allow us to separate a sufficiently clean and well-characterized bulge 
sample that we are able to detect a small population of bulge objects in the region 
of the color-magnitude diagram commonly occupied young objects and blue strgglers. 
However, variability measurements of these objects clearly establish that a fraction 
of them are blue stragglers. Out of the 42 objects found in this region of the color- 
magnitude diagram, we estimate that at least 18 are genuine BSS. We normalize the 
BSS population by our estimate of the number of horizontal branch stars in the bulge 
in order to compare the bulge to other stellar systems. The BSS fraction is clearly 
discrepant from that found in stellar clusters. The blue straggler population of dwarf 
spheroidals remains a subject of debate; some authors claim an anticorrelation between 
the normalised blue straggler fraction and integrated light. If this trend is real, then 
the bulge may extend it by three orders of magnitude in mass. Conversely, we find that 
the genuinely young (< 5Gy) population in the bulge, must be at most 3.4% under the 
most conservative scenario for the BSS population. 

Subject headings: Galaxy:bulge Galaxyrdisk - Galaxy:evolution 
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Introduction 



Blue Straggler Stars (hereafter BSS) are thought to be old, hydrogen-burning stars, made 
hotter and more luminous by accretion and now residing i n a region of the color-magnitude dia 



1993 



gram (CMP) normally occupied by much younger stars (e.g. lBond &: MacConnel]||l97ll : lLivio 
Strvkeiill993l). They have been d etected in examples of n early every stellar s ystem - open cluster s 



(e.g. lAhumada k. Lapasseti 120071 ). globular clusters (e.g . iPiotto et al 



stellar halo (Preston &: Sneden 



dwarf spher oidal (dSph) galaxies (IMomanv et al.l 120071: iMapelli et al.l 120071 ) and the Milky Way 



200C 



Carnev et al 



long been suspected (e.g. iMorrison k. Hardinglll993l ). their existence in the bulge as a population 



2004 



Lanzoni et al.l 120071 ) 



200 ll ). However, while BSS in the bulge have 



has never been proven due to strong contamination from disk stars that fill the same region in the 

cmdQ 

By the standards of modern space-baced astrometry, disk and bulge stars have large relative 
motions (~ 5 mas yr~^) and thus a reasonably pure sample of bulge stars can be isolated kinemati- 
cally. Prior use of the Hubble Space Telescope (hereafter HST) to perform this separation indicated 
that the population of bulge blue stragglers must be small (< few percent), but uncertainty in the 
kinematic contamination (and thus foreground disk objects mimicking BSS), combined with the 
absence of corroborati ng evidence like variability, prevented constraint on the BSS population size 
(iKuijken &: RichI |2002| ) . As a result, BSS have never unambiguously been detected in the Milky 
Way bulge. 

Beyond the intrinsic interest of BSS populations in a predominantly old stellar population, 

detection of BSS in the Milky Way bulge would be impor tant for galactic str uct ure formation. Most 

spira l galaxies show flattened "Boxy /Peanut" bulges (see Athanassoulall20o"5 and Bureau fc: Athanassoula 
20051 for a clear discussion) o r a small "Disky bulge" accompanied by ongoing star formation 
(jKormendv &: Kennicuttll2004l ) . Such a bar-driven arrangement of the sta rs evolves over a timescale 



of several Gy, and can be transient with multiple generations of bar (e.g. ICombesll2009ll. However 



chernical evidence suggests th e stars themselves mostly formed early and rapidly (e.g. iBallero et al 



20071 : iMcWilliam et al.ll2010l ). While the majority view now is that most stars in the Boxy /Peanut 
bulge of the Milky W ay (which we call here "the bulge") formed ~10 Gy ago (e.g. IZoccali et al 



2003 



FreemanI l2008l ) . the degree to which this star format i on was extended oy er time, is presently 



only somewhat weakly constrained (IKuijken &: Rich 



2002 



Zoccali et al.ll2008l ). Thus, better esti- 



mates of the young stellar population of the bulge, and thus its evolution, are a natural dividend 
of constraints on the BSS fraction. 

We report our use of a combination of proper motion and variability information to set this 
constraint. The observations are detailed in Section [2j Section [3] explains the selection of candidate 
BSS and the analysis used to estimate the size of the true BSS population. Finally in Sections |4]&: 



^The hypervelocity B-type star (HVS) HE 0437- 5439 is most likely a BSS from the inner milky way, although its 
BSS nature may result from the ejection-encounter ( Brown et al. 20ld ). 
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Owe discuss the BSS population in the context of stellar populations in general and the formation 
of the bulge of the Milky Way in particular. 



2. Observations 



All motions on which we report here were obtained from two-epoch imaging of the SWEEPS 
(Sagittarius Window Eclipsing Extrasolar Planet Search) field (Sahu et al. 2006) towards the 
bulge (at l,b = 2.65°, —1.25°) using the Wide Fi eld Channel of the Advan ced Camera for Surveys 
on HST, which has pixel scale 49.7 mas pix~^ ( van der Marel et al.l 120071 ). The sight-line passes 
through the bulge just beneath the galactic mid-plane, passing within ~350pc of the galactic center 



(but always exterior to the nuclear cluster, which shows ongoing star formation; iLaunhardt et al. 



20021) . Full details on the analysis of these extensive HST dat a sets c an be found in lClarkson et al 
(j2008l . hereafter Paper I), Sahu et al. (2006) and iGillilandl (j2008l ): we report here the aspects 
relevant to BSS detection and characterization. 



2.1. Astrometry and Photometry 



The measurements used are: (i) proper motions down to FSIAW = 24 in the SWEEPS field 
as measured from astrometry in the two epochs; (ii) absolute photometry in F606W &: FSliW, 
and (iii) time-series photometry of the stars from the continuous seven-day SWEEPS campaign, 
with gaps only due to Earth occultation. 

The techniques used to measure proper motions from the long (339 sec) and short (20 sec) 
exposures were detailed in Paper I. For objects saturated in the long exposures, motions and errors 
using the short exposur es are subst i tuted . The stellar time-series and absolute photometry were 
taken from the wor k of I Sahu et al.l ([20061), in which full details can be found. Difference image 
analysis techniques (IGillilandl l2008l . and references therein) were used to produce the time-series; 
the noise in the result is very close to Poisson for the majority of objects, and perhaps 30% higher 
than Poisson for saturated objects. Absolute magnitudes were produced from psf fitting to the 
image model against which variability was measured (Sahu et al. 2006). 



3. Sample selection 

The bulge population was extracted using cuts on motion and measurement quality, and the 
BSS candidates evaluated from this sample (Section 13. Ij) . along with contamination by disk and 
halo stars that fall into the sample (Section 13. 2p . Variability information was used to set limits on 
the size of the young (< 5 Gy) population among bulge BSS candidates (Section l3.3p . corrected for 
periodicity detection incompleteness (Section l3.4p . This allowed the size of the true BSS population 
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in the bulge to be bounded (Section l3.5p . 

3.1. Separation of bulge from field 

The streaming of disk stars in front of the bulge in galactic longitude was used to isolate the 
bulge population. Proper motions were measurable for 188,367 of the 246,793 stars with photometry 
in the field of view. Stars were selected for further study with proper motion error ai < 0.3 mas 
yr~^ and < 0.45 mas yr~^, and with crowding parameter q < 0.1 (see Paper I for definition of 
the crowding parameter). This latter limit is slightly more restrictive than used in Paper I, and 
results in a well-measured astrometric sample of 57,384 objects down to FSIAW ~ 23. Discarding 
objects saturated in the short exposures then leaves 57,265 well- measured, unsaturated objects. 
To set the proper motion cut for bulge stars, the distribution of latitudinal proper motion fii was 
measured for stars above the main sequence turn-ofi' (hereafter MSTO), where the disk+BSS and 
bulge populations trace separate loci in the CMD. Using a proper motion cut /i; < —2.0 mas yr^^, 
rejects all but 0.19% of disk stars, while admitting 26% of bulge objects. With error cut at 0.3 mas 
yr~^ and ^/ < —2.0 mas yr^^, galactic-longitudinal motions must be detected to at least 6a for 
inclusion in the cleaned bulge catalogue. Applied to the full error-cleaned sample with astrometry, 
the motion cut produced a population of 12,762 likely-bulge objects. 

The BSS selection region was set with the constraints that (i). the bright end should not 
impinge on the region in which blue horizontal branch objects are expected, and (ii). the faint and 
red ends should not greatly overlap with the main population of the bulge near the MSTO (see 
Figure 1 for the selection region adopted in the F8UW, F606W - F8UW CMD). Within this 
selection-region, 42 objects were found in the kinematically-selected bulge sample. 

3.2. Kinematic and photometric contaminants 

The number of disk objects amongst our BSS candidate sample was estimated directly from 
the distribution of longitudinal proper motion fn. Beginning with the 57,265 objects that are 
astrometrically well-measured, the longitudinal proper motion distribution is assessed within the 
post-MSTO bulge sample as shown in Figure [2] (top panel). The CMD selection-region for this 
post-MSTO bulge sample was constructed by inverting the BSS selection box in color about its 
reddest point, while keeping the same magnitude range as the BSS sample. We denote the set of 
longitudinal proper motions of these post-MSTO bulge objects as fii^buige- It is clear from Figure 2 
that the contribution of any disk objects to this population is very small. This population was fit 
with a single gaussian component. 

The longitudinal proper motion distribution within the BSS selection region we denote by 
f^i,BSSj was fit with a model consisting of two gaussian components, one (fJ-i^i) representing 
disk objects, one (/U/^2) representing bulge objects (Figure O middle panel). The centroid and 
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width of were frozen to the values fit from Hi^buige', because the BSS and post-MSTO bulge 
samples cover an identical magnitude range, to first order they should suffer similar instrumental 
and incompleteness effects. The expected number of disk objects within the BSS selection-box 
(denoted A''^) was evaluated from the parameters of the best fit to hi^bss with these constraints 
(the best-fit to Hi^bss shows xl = 24.2/26). 

We performed a Monte-Carlo simulation where a large number of proper motion datasets 
drawn from the best-fit model were generated, and the distribution assessed of recovered values 
of both the number of disk objects in the BSS region and the number that would pass our 
kinematic cut for bulge objects (-/V^,<-2.o)- The results are shown on the bottom panel of Figure [2j 
The recovered -/V^,<-2.o is 0.7 it 0.73 objects, but the distribution is rather asymmetric - while the 
most common recovered suggests fewer than one contaminant, in some cases up to four 

contaminants are measured. We therefore adopt 0-4 as the ranges on the contribution from disk 
objects to our kinematically cleaned bulge sample. 

Of the ~ 6 estimated halo contaminants (Paper I), of order 0-2 objects would fall into the BSS 
selection region. So, of 42 objects, we expect — 6 objects to have passed our kinematic cuts that 
do not belong to the bulge. 

We investigated by simulation the contamination in our BSS region due to the main bulge 
population near the MSTO, which might be expected to put o utliers in the BSS r egion. Synthetic 



bulge populations were constructed with age range 7-13 Gv (IZoccali et al 



2OO3I ) and metallicity 



distribution approximating the bulge ( Santiago et al.ll2006l : IZoccali et al.ll2008l ). Age and metallicity 
samples were conv erted to predicted instrumental magnitudes using the isochrones described in 
Brown et al.l (j2005l ). The resulting population was then perturbed by the bulge distance distribution 
and with a bina ry population est imated using a variety of binary mass-ratio distributions and 
binary fractions (jSoderhjelml 120071 ) to produce simulated bulge-selected populations in the CMD. 
The resulting number of objects from the main bulge population that happen to fall in the BSS 
region in the CMD, is 5 it 2, over 10^ trials. 

We therefore estimated the total contaminant population at (5-13)/42 objects, leaving (29- 
37)/42 possible bulge BSS. 



3.3. Photometric Variability 



BSS may form by a variety of mec hanisms (e.g. lAbt 



1985 



Livid I1993I ). including mass trans- 



fer in a binary system (IMcCreal Il964l'). This can yield BSS that are presently in bin aries (e.g. 



Andronov et al.ll2006 ; 



C 



2001: Mathieu fc Celler 



len fc Hanll2009l ') with observable radial- velocity periodicities (e.g. lCarney et al. 
2OO9I ). For periods <10 days, photometric variations caused by tidal defor- 



mation of one or both mem bers may be observed, in the extreme producing a W UMa-type contact 
binary (lMateolll993l . Il996h . 
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The variability information from the 2004 epoch (Section [2]) allows the tendency for W UMa 
variability among the BSS candidates to be examined. In a manner similar to Section I3.H a 
population composed of mostly disk objects was isolated, this time using longitudinal proper motion 
m > +3.0 mas yr~^ and again requiring 6cr motion detection. This test sample indicates the fraction 
of W UMa-like variables among a young population observed at the same signal to noise range as 
the bulge objects we wish to probe. Comparison of the W UMa occurrence rate between the 
mostly-disk and mostly-bulge samples then provides an estimate of the fraction of genuine BSS 
among the BSS candidates. 

To avoid confusion with pulsators, systems were not counted among the W UMa if the dominant 
period detected P was shorter than 0.5d , or if significant phase or frequency variation was present 
along the observation interval (iGillilandl 120081 ) . Note that we do not eliminat e possible pu lsators 
from our BSS candidate table, as BSS can also show photometric pulsations (jMatedll993l ). Each 
variable lightcurve was visually examined to determine if its morph ology resemb l es that of W 



Uma objects or close eclipsing binaries found in other systems (e.g. iMated Il993l : lAlbrow et al 



2001 



). W UMa's were much more populous amongst the bulge-selected population: of the 42 
bulge-selected objects, four show W UMa-type variability (one additional object is apparently a 
long-period pulsator; Figure [3|). Of this four, one object only has a single eclipse event recorded and 
awaits confirmation. Conversely, none of the 81 disk-selected objects show W UMa photometric 
variability. 



3.4. Periodicity completeness correction 

The completeness of the Lomb-Scargle (hereafter LS) variability search to W UMa variables 
was estimated through simulation. Real lightcurves with no significant periodicity detect ion (peak 



LS p ower <5) were injected with W UMa-type lightcurves (following the description of iRucinski 



19931 ) under a log-uniform period distribution. Synthetic W UMa systems were denoted "detected" 
if the peak in the Lomb-Scargle periodogram was within 10% of the input period, the signal was 
present in the lightcurves in both filters, and showed LS peak power greater than 13 (the 99% 
significance level, defined as the peak LS power exceeding that returned by 99% of trials with no 
periodicity injected). The process was repeated for input variation amplitudes Amag/mag = 0.1, 
0.05, 0.02, 0.01, 0.005, 0.002, 0.001 and period bins 0.4, 0.8, 1.7, 3.4 and 7.0d (Figure [3]). The 
dependence of the detection completeness on injected period is rather weak but the dependence on 
amplitude is clear; at least 88% of systems with 2% fractional amplitude variation were detected, 
and 70% at the 0.5% amplitude level. 

To i nap signal amplitude to sys tem inclinations, synthetic W UMa lightcurves were simulated 
using the Wilson Sz Devinnev ( 1971 ) code as implemented in the Nightfall packag^ indicating that 



^http: / /www.hs.uni-hamburg.de /DE /Ins/Per /Wichmann/Nightfall.html 
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at most -^70% of W UMa type systems would be detectable down to binary inclination ~45°, so 
we are sensitive to perhaps 35% of the true W UMa population. This implies that the 3-4 detected 
W UMa trace an underlying population of 9-11 systems with W UMa-like system parameters, so 
at least 9-11 of the 42 BSS candidate objects have W UMa system parameters, and are therefore 
not young. 

If, as our measurements suggest, 9-11 of the 42 BSS candidates are present-day W UMa 
objects, we can ask what fraction of the BSS candidates are in binaries of any orbital period. The 
period distribution for bulge BSS presently in binaries is unknown; we use observations of the solar 
neighbourhood and open clusters as proxies for the bulge population. 



Carney et al.l (120011 ) report spectroscopic periods for 6/10 high-proper motion blue metal-poor 
(hereafter BMP) stars in the solar neighborhood, isolating BSS candidates by virtue of de-reddene d 
colors bluer than the MSTO of globular clusters of comparable metallicity (ICarnev et al.lll994l) 



That study also reports a re-analysis of the radial velocities of the BMP sample of lPreston &: Sneden 



(|2000l ). finding ten spectroscopic BMP orbits for a grand total of sixteen spectroscopic orbits 
among 6 BSS candidates and 10 further BMP that may be BSS. This yields 2/16 objects with 
logiQ{Porb) < 1 and 14/16 with 1.0 < logio(P) < 3. 2 (with one possible very long-period binary 
not included in this sample). Thus, in the solar neighborhood perhaps 1/8 of the total binary BSS 
population resides in short-period binaries, with log^ol-^orb) < 1-0. 



Turning to the open clusters, iMathieu &: Gelled (|2009l ) find a high binary fraction and broadly 
similar period distribution amongst BSS candidates in the ~ 7 Gy-old open cluster NGC 188. There, 
16/21 BSS candidates are in binaries. Of these sixteen objects, two show log^gl-P) < 1-0 while the 
remaining 14 a ll show 1.0 < logio(-P) < 3.5. For the ~ 5 Gy-old open cluster NGC 2682 (M67), 



LathamI (j2007l ) r eports two BSS in binari es with log]^o(-for6 < 1) and five with 1.0 < logiQ{Porb) ^ 
3.1 (Figure 3 of IPerets fc Fabrvckvl 120091 ) . Thus, in open clusters the short-period binaries may 



< 



make up between about 1/8 and 2/7 of the total population of BSS in binaries. 



The period distribution for main-sequence binaries of several types is collated in lKroupal (jl995l ). 
For < 3.5 the numbe r of binaries increases with P^^^ at least as steeply as a log-uniform 
distribution (jKroupal Il995l ). If BSS in binaries in the bulge follow a similar distribution, then 
binary BSS with logi(^{P) < 1 would make up at most about 1/3 of the total population with 
logio(P) < 3.0. 

Taken together, these estimates suggest that, if all the W UMa are indeed BSS presently in 
binaries, then BSS in binaries could by themselves account for all the 29-37 non-contaminant objects 
among our candidate BSS. 
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3.5. Bounds on the true BSS population in the bulge 



As none of the kinematically-selected disk objects show W UMa characteristics, the W UMa 
stars in the BSS region are unlikely to be disk stars. Furthermore, W UMa variability is expected 
to be a natural evolutionary state of some of the blue stragglers. Thus, the detection of W UMa 
variability in the BSS region of the bulge CMD certifies these stars as true bulge blue stragglers. 
Hence, we conclude that the 9 — 11 W UMa objects are in fact bulge BSS, from a possible bulge 
BSS population of 29-37 (Section 13. 2p . This leaves 18-28 bulge objects whose nature remains 
undetermined. It is unlikely that W UMa variables make up the entire BSS population in the 
bulge. Based on arguments in Section 13.41 we consider the BSS population with binary periods 
log^Q(P) < 1 to comprise at most half of the population of BSS with any orbital period. We 
therefore consider two extreme scenarios: 1. Optimistic: all of the remainder are bulge BSS thus 
there are no genuinely young bulge objects in our sample, leaving 29 — 37 BSS, and 2. Conservative: 
our W UMa objects make up half of the BSS among our sample, suggesting 18-22 genuine BSS and 
7 — 19 possibly young bulge objects. 

In addition, a significant fraction of the bulge BSS population may exi st with no binary com 



panion or m 



2006 



very 



wide binaries, as indicated by population studies (e.g. Andronov et al. 20od; 



Chen &: Hani |2009|) and observations of clusters and the solar neighbourhood (e.g. 



Ferraro et al, 



Carnev et al.l 12001 



Therefore the true BSS population in the bulge probably lies between 
the "conservative" and "optimistic" scenarios outlined above. 



Discussion 



4.1. Normalization of Nbss 



As probes of binary evolution, BSS indirectly probe the star formation history of stellar popu- 
lations in their own right. To compare the bulge BSS population with that of other stellar systems, 
we normalize by the number of horizontal branch objects in the system; Sbss = ^BSS /^HB- We 
adopt Nhb = Nbhb + Nrhb + Nrc, where the three terms denote the Blue and Red horizon- 
tal branch stars, and the Red Clump stars respectively, and thus cover the full metallicity range 
of post-MS core H e -burning objects. This is the convention used for the dwarf spheroidals (e.g. 
Carrera et al.ll2002l : iMapelli et al.ll2007l : iMomanv et al.ll2007l . although som e authors choose not to 
include Nrq explicitly if none are fourid in their sample) ; open clusters (jde Marchi et al.l |2006| ) 
and the globular clusters (jPiotto et al.ll2004l . though here the definition of Nhb is not explicitly 
stated). The exception is the stellar halo, where Nhb = Nbhb is adopted because only BHB can 



be di stinguished kinematically and photometrically from the local disk population (jPreston et al 
199lh . 



The kinematic cuts leave too sparse a population to accurately disentangle the RC from the 
underlying Red Giant Branch (RGB). We therefore estimate Njic from the F606W histogram along 
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the RGB only after selection by astrometric error and crowding parameter, and scale the result to 
the bulge population passing the full set of cuts (Section 13. ip . Figure H] shows the selection; the 
marginally-detected RGB Bump (hereafter RGBB) is masked out as well as the RC before fitting the 
continuum due to the RGB; Gaussian components fit to the RC and RGBB then yield the number 
of objects in each population over the RGB continuum. Bootstrap monte carlo sampling before 
producing the histogram is used to estimate the fitting uncertainty for a given binning scheme. We 
find Nrc = 180 ±53 while the RGBB is detected at more marginal significance; Nrgbb = 41 ±2311 
Our estimate of Nrc = 180±53 scales to the 26% of bulge objects surviving the full kinematic cuts, 
as 47 ± 14 objects. Even in the short exposures, roughly 15% of this population would saturate 
and thus not enter the bulge sample (Figured]), so that we would expect 40 ± 12 RC objects to fall 
within the bulge sample. This figure includes the RHB objects. 

The contribution from BHB is difficult to estimate for the bulge, but is clearly rather small - 
it vanishes against the much larger foreground disk population without kinematic selection. Unlike 
the RC±RHB selection region, saturation of the astrometric exposures is not an issue for the 
BHB selection region we adopt (Figure [1]). Altogether, four objects are found in the BHB ± 
RHB selection region; of order one object may be a kinematic contaminant. We therefore consider 
(3 — 4 ± IcTpoisson) = 2 — 6 as the limits on Nbhb- 

Thus we estimate the total Nhb =30-58 objects among the same sample from which the BSS 
were estimated. Our estimate of 18 < Nbss ^ 37, then leads to 0.31 < {Sbss = ^BSS /^hb) < 
1.23. 



4.2. Comparison to other systems 



For the halo, iPreston k. SnedenI (|2000l ) estimate Nbss/^bhb = 5, though their selection 
criteria appear to admit a population that is a factor 2 lower; we adopt 2.5 < Sbss ^ 5 for the 
halo. The upper limit Nbhb = 6 for the bulge indicates the lower limit Nbss/Nbhb > 3.0. Thus 
our BSS fraction is marginally consistent with that of the halo, but as yet the statistics are too 
poor to support further interpretation. 

In principle it would be of interest to compare the normalised blue straggler population in the 
bulge with those from the dSph and the clusters. In partic ular, the collect i on of Sbss vs integrated 
light My (a proxy for the mass) for clusters and dSph of iMomanv et al.l (|2007l ) suggests differing 
trends for clusters and spheroidals, against which the bulge would make an interesting comparison. 



Nataf et al.l (|201ir ) recently found an anomalously small RGBB population compared to the RC, with 



Nrobb /Niic = {^'2.7 ±2)% among a large number of sight-lines across the bulge. We do not confirm this finding (com- 
pare our Figure 4 with Figure f of lNataf et aLll2011f ). although with Nrgbb /NRc = i23 ± 15)% our statistics are too 
poor to falsi fy their claim of He-enrichment. In addition, s ampling as we do towards the center of the "X" structure 
in the bulge ( Nataf et al. 2O10l : Mc William fc Zoccalil 2010|). an intri nsic difference in Nrgbb/Nrc may be expected 
between our sight line and the range reported bv lNataf et al.l (l201l|). 



-lo- 



in practice, however, the true variation of blue straggler population with mass, is still far 
from settled observationally; we use the case of the dSph objects Draco and Ursa Minor as a 
case study. For these objects, different authors reach strikingly different conclusions using data 
from the sar ne facility (the Wide Field Camera on the 2.5m Isaac Ne wton Telescope). In their 
compilation, Momanv et al.l (120071 1 quote Sbss ~ 1. 23,1.35 for D r aco (Carreraet al.ll2002l ) and 



Ursa Minor (jAparicio et al 



2001 



respectively, while iMapelli et al.l (120071 ) obtain instead Sbss ~ 
0.25, 0.3 respectively for the same systems. (Both studies conclude that Draco and Ursa Minor are 
true fossil systems, with no extensive present-day star formation activity.) 



The two studies differ in the choice of selection region for BSS on the C MP; that of lMapelli et al 



(j2007l ) is about a factor two smaller than that of iMomanv et al.l (|2007l . and refs therein), and is 
separated from the MSTO by a larger distance in the CMD. Thus it is unsurprising that the more 
conservative BSS region produces a lower number of BSS candidates; our problem is choosing the 
study most appropriate f or comparison t o our o wn work. The BSS selection region we have adopted 
is most similar to that of lMomany et al.l (|2007l . compare our Figure 1 with Figure 1 of that paper), 
and is consist ent with the sel e ction regions used for the cluster studies cited therein. Our belief is 
therefore that I Momanv et al.l (|2007l ) is the most appropriate for comparison to our own estimate. 



In addition to photometric selection effects, the sample of objects under consideration includes 
systems of differing turn-off mass, which means the BSS populations uncovered in these systems 
may have quite different formation mechanisms from each other. Furthermore, our sample of 
BSS candidates is different from those of the external systems since our candidates were found in a 
pencil-beam survey through a narrow part of the bulge. A better interpretation requires population 
synthesis modeling, which is beyond the scope of the present in vestigation. Wi t h the se caveats in 
mind, we present the BSS fraction of the bulge compared to the lMomanv et al.l (|2007l ) selection in 
the hope that it will provoke just such an investigation. 



We use COBE photometry to estimate My ~ 20.4 for the bulge (jDwek et al.lll995l ). When 
placed on the Sbss I M y diagram, the bulge appears to be c onsistent with the tren d pointed out by 



Momany et al 



20071 ) . If the {My, Sbss) trend reported in Momany et al. (j2007l ) is indeed borne 
out by further observation, then the bulge may extend this trend by over seven magnitudes in 
My (or, assuming constant M/L ratio, nearly three orders of magnitude in mass; Figure [5]). What 
this suggests about BSS evolution in the bulge as compared to the dwa rf spheroidals, is unclear 
at present, and we defer interpretation until the iMomanv et al.l (|2007l ) trend has been further 
established or falsified. It should also be borne in mind that our study samples the bulge quite 
differently (narrow and deep) to the dSph galaxies (wide-area but relatively shallow); however as 
ours is the first study to provide a measurement of Sbss for the bulge, we show its location on the 
Sbss-^v diagram to stimulate further investigation. 



We can say that the bulge is discrepant from the sequence suggested by the globular clusters; 
for it to lie on the cluster trend we would have to have observed no BSS in our sample at all. 
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4.3. Young-bulge population 

For the purposes of this report we define "young objects" as stars already on the main sequence, 
with ages < 5 Gy. Our photometry covers main sequence objects up to about two solar masses, so 
the main sequence lifetime of all objects in our sample is long compared to the time spent on the 
Hyashi track. 

The simplest estimate assumes that the fractional contribution by young objects is the same 
across all luminosity bins. In this case the young-bulge population is estimated from the number of 
non-BSS bulge objects in the BSS region of the CMD (0-19; Section [33]) and the total number of 
bulge objects in the same magnitude range and with the same kinematic and error selection (346 
objects). This yields a young-bulge fraction in this luminosity strip, of 0-5.5%. 

A more complete estimate extrapolates the young population within the BSS selection-box 
to the entire bulge sample for which we have astrometry. This extrapolation is complicated by 
differing incompleteness to astrometry above and below the MSTO (Section [2]), and by our lack of 
knowledge of the luminosity function of young bulge objects. 

To estimate this scaling, we make two assumptions: (i). that the young bulge population and 
disk population both follow the same Salpeter mass function (i.e. dN/dm oc m~'^'^^) and the same 
mass- magnitude relationship, and (ii). that the astrometric selection effects for disk and young 
bulge objects vary with instrumental magnitude in the same way. The second assumption allows 
us to account for astrometric incompleteness in the young-bulge objects, while the first allows us to 
correct the estimate for the fact that the disk and bulge populations are seen at different distance 
moduli and thus sample different parts of the mass function. 

Denoting the number of objects within the BSS selection region as Nbss and the number 
from the base of the BSS region in the CMD to the faint limit of our kinematically-cleaned sample 
(F814VF = 23) as N faint we expect 



/ N faint \ _ f C faint \ f nfgjnt \ 

V Nbss J disk V cbss J disk V ^bss J disk 



f N faint \ _ f C faint \ f ^ faint \ ^^^^ 

V Nbss ) young bulge V ^BSS J young bulge V ^BSS ) young bulge 



where the quantities cbss and ubss refer to the astrometric completeness and the integral of 
the mass function respectively within the magnitude limits of the BSS selection box. We do 
not distinguish between the initial and present-day mass function for the disk and young bulge 
populations. 



Under our assumption (ii) above that (l^^) = ^ estimates of f^^^) 

\ BSS J disk V SSS / young bulge V SSS / disk 

and f^^) allow us to solve for f^^) . The first two quanti- 

V SSS / yQung bulge V SS'S / disk V SSS / young bulge 

ties may be estimated from the mass function under assumption (i) above, while the third can be 
estimated from our kinematic dataset. 
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We use the zero-age, solar metallicity disk isochrone of lSahu et alj (|2006l ) to convert magnitude 
limits to mass limits for both the disk and young-bulge populations. Investigation of a wider range 
of metallicities for young objects is beyond the scope of this report. In reality, the disk stars are 
distributed over a range of distances along the line of sight, although as the largest contribution of 
disk stars comes from the Sagittarius spiral arm we restrict ourselves here to a single distance for 
the disk population. The isochrone allows us to convert magnitude limits in the "BSS" and "faint" 
samples to mass limits for the disk population in order to estimate I ) = fMFdisk- To 

estimate this quantity for the young bulge population, the difference in distance moduli for bulge 
and disk are used to estimate the corresponding mass limits for the young-bulge sample. Table 2 
shows the mass limits adopted and the correction factors adopted from this analysis. 



All that remains is to estimate 



from observation. We follow a similar method to 



^BSS / disk 

the fitting of the size of the bulge and disk population within the BSS region as described in Section 
l3.2l above: the longitudinal proper motion distribution fj^i is estimated for a series of magnitude strips 
beneath the MSTO and the number of disk objects in each magnitude bin is estimated from the 
component of the fit that corresponds to the disk population. The results are in Figure [6l the 
number of disk objects in the BSS region is N^i bss = 354 it 23 while the number below the BSS 

N faint \ 
. ^BSS J disk 

' N faint \ _ ^o/i Q _i_ o /i\ /'5.20^ 



24.8 ±2.4. This then 



region is Ndjaint = 8781 ±617 objects, resulting in the scaling 

(24.8 ± 2.4) X fl^) = 22.0 ± 2.1. We remind the reader that we are 



means 



Nbss 



young bulge 



^5.84; 



concerned in this section entirely with the young bulge objects. 

This means that under our "conservative" scenario, the 7-19 young bulge objects suggested 
by our dataset within the BSS region would scale to approximately 163-437 young bulge objects 
across the entire magnitude range. As a fraction of the 12,762 kinematically-selected bulge objects, 
this young population would make up as much as 1.3% — 3.4% of the bulge population, under our 
"conservative" scenario for the number of true BSS in the bulge. 

The small young-bulge population is difficult to reconcile with the recent discovery of a metal- 
licity transition (high to low movi ng outward from t he galactic center) that has been reported at 
~0.6-lkpc beneath the mid-plane (jZoccali et al.ll2008l ). If this transition really does indicate a more 
metal-rich population interior to 0.6kpc, then the lack of young bulge objects along our sight line 
would suggest t hat this metal- r ich in terior population is not a younger population. Corroborating 
this suggestion, lOrtolani et al.l (jl995l ) considered the age of the bulge in Baade's Window and found 
it similar to that of the globular clusters. 

In addition, a growing number of main sequence objects in the bulge are becoming spectro- 
scopically accessible by microlensing; at present fifte en bulge dwarfs hav e been spectroscopically 
studied in a box 10° x 8° about the galactic center (|Bensbv et al. Il20ld ). Among these objects, 
three are spectroscopically dated to < 5 Gy, far more than the < 3% we would predict. If we have 
somehow vastly underestimated our kinematic contamination from the disk, the discrepancy is only 
sharpened because the true sample of young bulge objects would then be reduced still further. At 
least part of this discrepancy may be due to the fact that most of the lensed sources are not in the 
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inner bulge region that we probe here. In addition, star formation rates within the b ulge probably 



do v ary strongly over the large spatial scales probed by microlensing studies (e.g. iDavies et al 



20091 ). We note that all three young objects in the Bensby et al. (2010) sample (their objects 310, 



311 and 099) are degrees away from our field on the plane of the sky, at negative galactic longitudes. 
We await further results on the spatial distribution of metallicity and age within the bulge, with 
great interest. 



5. Conclusions 

HST observations of a low-reddening window through the bulge have yielded the first detection 
of Blue Straggler Stars (BSS) in the bulge of the Milky Way. By combining kinematic discrimination 
of bulge/disk objects with variability information afforded by an HST dataset of unprecedented 
length (from space and for the bulge), we find that: 

• Of the 42 BSS candidates identified with the bulge, between 18-37 are true BSS. We estimate 
for the horizontal branch Nhb = 30-58 with the same kinematic selection, suggesting 0.31 < 
Nbss/Nhb < 1.23. 



If t he trend in normalized BSS fraction against integrated light suggested by iMomanv et al 



(|2007l ) is real, the bulge may extend this trend by nearly three orders of magnitude in mass. 

• Normalized appropriately, the BSS population in the bulge is marginally consistent with that in 
the inner halo but the allowed range for the bulge is very broad. 

• The truly-young (< 5Gy) stellar population of the bulge is at most 3.4%, but could be as low as 
along our sight-line. 



If the recently-discovered metallicity transition ~ 0.6 — Ikpc beneath the mid-plane (jZoccali et al 



20081 ) does indeed indicate a population transition, then the more metal-rich population does not 



have a significant component of age < 5 Gy. 

In the old metal rich Galactic bulge, we may therefore conclude that blue stragglers are a 
small component of the stellar population (18-37 of 12,762 kinematically selected bulge objects), 
and would not significantly affect the integrated light of similar unresolved populations. 
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Fig. 1. — Color-Magnitude Diagram (CMD) of the kinematically selected extreme-bulge population. 
Isochrones are shown for a young foreground population (blue dash-dot), old bulge population 
(dashed) and an intermediate-age, slightly metal-poor bulge population (solid). Diamonds show 
the median bulge population. Ages are indicated, as are the adopted [Fe/H] values. Selection 
regions for BSS and horizontal branch objects are indicated. Saturation limits in the short and 
long exposures are marked, and the five BSS candidates exhibiting close binary variable (W UMa) 
lightcurves are indicated. See Section [3?T1 
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Fig. 2. — Estimation of the kinematic contamination in the BSS CMD-selection region (Figure [T]), 
due to disk objects. Top: single-component gaussian fit to the distribution of longitudinal proper 
motion for post-MSTO bulge objects (nifiuige), which we use to calibrate the bulge component 
in the BSS region. Middle: two-component gaussian model fit to the longitudinal proper motion 
distribution within the BSS region of the CMD (m^BSs)- Bottom-Left: distribution of recovered 
number of disk objects in the BSS region under 5,000 Monte Carlo trials. Bottom Right: Distribu- 
tion of recovered number of disk objects that would pass the kinematic cut /i/ < —2.0 mas yr~^. 
See Section [321 
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Fig. 3. — Folded lightcurves for the five objects sliowing sinusoidal or eclipse-type variability with 
a period longer than 0.5 days. The object at the top-right is probably a long-period pulsator. The 
objects at middle-left and bottom-left were detected at half the likely true period, as suggested by 
differing eclipse depths half a period apart. Bottom Right: Completeness to W UMa variability 
from the HST dataset. Left: fraction of injected W UMa lightcurves recovered as a function of 
magnitude and amplitude. Injected variation amplitudes AV/V are labeled next to each set of 
curves. For each set of curves, from top to bottom, symbols represent periods 0.4, 0.8, 1.7, 3.4 and 
7.0d respectively; in practice, all the limits for a given amplitude overlap except for the 7d set. 
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Fig. 4. — Estimation of the red clump population Nrc for the bulge. We use the sample selected for 
astrometric error and crowding, and use the population within the bulge RGB rather than cutting 
by proper motion to afford sufficient statistics to constrain Nj^c accurately. Diamonds show bins 
selected for the fit to the underlying RGB population (dashed line), subtracted before estimating 
Nrc- The fits to the RC population (solid line) and RGB bump (dotted line) are indicated, as 
are the population sizes in the RC and RGB bump (the latter is marginally detected). Triangles 
indicate astrometric saturation limits for the RC at the blue and red ends of the selection box 
(Figure 1). See Section HTT] 
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Fig. 5. — Comparison of the normalized bulge blue straggler fraction Sbss (limits: black diamonds) 
with other stellar populations, as a function of integrated light My. These are: dwarf spheroidal 
galaxies without prominent current star formation (asterisks; each is labeled), Milky Way globular 
clusters (boxes), the Milky Way globular clusters NGC 6717 and NGC 6838, which may have sig- 
nificant field-star contamination (triangles) and the Milky Way open clusters (crosses). Literature 
points are from the compilation of Momany et al. (2007 and refs therein) except for the bulge (this 
work), and the redeterminations of the BSS fraction in Draco and Ursa Major by Mapelli et al. 
(2007; indicated with arrows as they overlap the cluster sequence). See Section [4.21 
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Fig. 6. — Estimating the scaling of the number of young objects within the BSS selection region 
(Figure 1) to the full kinematic sample, in the presence of datasets of differing depth. For each 
magnitude strip, the marginal fii distribution is fit by a two-component gaussian, identified with 
the disk and bulge. The fractional contribution from the disk is then identified from the fits, and 
the errors estimated by Monte Carlo (Section 13. 2p . Top: disk fraction among the magnitude strips 
below the BSS CMD region; the brightest bin contains the MSTO itself. Bottom: histogram of the 
number of disk objects detected in the sample, this time with the number of disk objects detected 
in the BSS region included (brightest bin). See Section [4.31 
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Table 1. Table of Bulge Blue Straggler Candidates. 



f814W f606W RA (J2000.0) Dec (J2000.0) P(d) 

16.397 17.458 17^ 58™ 55.441* -29° 11' 1.110" 

16.448 17.426 17^ 59™ 1.516" -29° 10' 49.401" 1.152 

16.537 17.573 17^ 58™ 58.208*' -29° 12' 1.974" 

16.544 17.524 17^ 59™ 7.983* -29° 11' 8.239" 

16.705 17.811 17^ 59™ 5.742* -29° 12' 3.695" 

16.822 17.941 17^ 59™ 2.823* -29° 10' 18.616" 

16.921 17.885 17^ 59™ 6.525* -29° 13' 25.552" 

16.942 17.955 17^ 58™ 59.127* -29° 12' 56.727" 

16.975 17.908 17^ 59™ 1.232* -29° 11' 14.057" 

17.121 18.157 17^ 59™ 1.975* -29° 11' 44.373" 

17.175 18.111 17^ 59™ 7.283* -29° 12' 56.907" 

17.207 18.105 17^ 58™ 57.713* -29° 12' 0.897" 

17.237 18.275 17^ 59™ 4.579* -29° 11' 1.411" 8.578 

17.330 18.316 17^ 59™ 6.396* -29° 10' 23.907" 

17.345 18.378 17^ 59™ 2.214* -29° 10' 45.657" 

17.454 18.446 17^ 59™ 4.413* -29° 10' 28.276" 

17.500 18.503 17^ 58™ 55.138* -29° 13' 4.462" 

17.512 18.598 17^ 59"^ 7.624^ -29° 10' 32.762" 

17.515 18.349 17^ 58"* 58.896* -29° 13' 15.047" 

17.537 18.507 17^ 59"* 1.420^ -29° 12' 21.630" 1.152 

17.552 18.507 17^ 58"* 53.865* -29° 10' 32.655" 

17.554 18.496 17^ 58"* 54.981* -29° 11' 18.907" 

17.557 18.552 17^ 59"* 3.068^ -29° 12' 52.866" 

17.559 18.444 17^ 59"* 5.433* -29° 12' 10.876" 0.681 

17.559 18.637 17^ 58™ 56.862* -29° 11' 8.062" 4.084 

17.564 18.571 17^ 58™ 53.759* -29° 13' 7.067" 

17.586 18.609 17^ 59™ 8.024* -29° 11' 35.341" 

17.604 18.527 17^ 59™ 6.340* -29° 11' 16.109" 

17.646 18.647 17^ 59™ 4.582* -29° 11' 46.148" 

17.680 18.659 17^ 58™ 55.199* -29° 11' 22.063" 

17.744 18.694 17^ 58™ 58.157* -29° 10' 54.293" 

17.783 18.833 17^ 58™ 59.460* -29° 10' 37.722" 

17.788 18.697 17^* 59"* 2.043* -29° 10' 19.798" 
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Table 1 — Continued 



f814W f606W RA (J2000.0) Dec (J2000.0) P(d) 



17.799 18.839 

17.805 18.764 

17.853 18.794 

17.853 18.868 

17.865 18.906 

17.877 18.915 

17.915 18.933 

18.106 19.110 

18.131 19.061 



17f 59™ 0.994* 
17^ 59™ 7.892* 
17^ 59™ 6.990* 
17^ 58™ 56.398* 
17^ 59™ 3.512* 
17^ 59™ 2.502* 
17^ 58™ 53.989* 
17^ 59™ 6.057* 
17'* 58™ 55.690* 



-29° 12' 2.624" 

-29° 10' 50.466' 

-29° 10' 28.606' 

-29° 12' 16.750' 

-29° 13' 22.640' 

-29° 13' 26.309' 

-29° 10' 26.682' 

-29° 11' 18.876' 

-29° 12' 34.222' 
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Mass limits - BSS region 
(16.2 < F8UW < 18.2) 


Mass limits - "Faint" region 
(18.2 < F8UW < 23.0) 


/ faint \ r 


Disk 
Bulge 


1.12 < # < 1.68 

— Mq 

1.41 < IT- < 2.11 

— Adp, 


0.43 < IT- < 1-12 

— Mq 

0.53 < 7^ < 1-41 


fMF,disk = 5.84 
fMF,bulge = 5.20 



Table 2: Mass-function correction used when scaling the number of objects within the BSS selection 
region to the full magnitude range (16.2 < FSIAW < 23.0). See Section for details. 



